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Abstract 

We propose a di-interstitial model for the P6 center commonly observed in ion 
implanted silicon. The di-interstitial structure and transition paths between 
different defect orientations can explain the thermally activated transition of 
the P6 center from low-temperature Cih to room-temperature D2d symmetry. 
The activation energy for the defect reorientation determined by ab initio cal- 
culations is 0.5 eV in agreement with the experiment. Our di-interstitial model 
establishes a link between point defects and extended defects, di-interstitials 

providing the nuclei for the growth. 

PACS numbers: 61.72.-y, 61.72.Cc, 61.72.Ji, 71.55.-i 
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Transient enhanced diffusion (TED) in boron-implanted silicon is the limiting factor in 
controlling dopant profiles for submicron Si-based devices. Interstitial defects in bulk Si 
generated during implantation have been identified as the sources for boron TED 
A class of macroscopical interstitial defects, namely {311} defects, was suggested to emit 
interstitials that can contribute to the enhancement of boron diffusion under typical implan- 
tation conditions ||l],@]. Decreasing the ion implantation energy can suppress the formation 
of macroscopic {311} defects, thus reducing boron TED. However, boron TED exists even 
at very low implantation energy (below 10 keV) in samples with no visible {311} defects 
0. This implies that microscopic interstitial defects, e.g., interstitial clusters, contribute to 
boron TED. 

The activation energy of boron TED, a measure of the energy required to dissociate 
interstitials from the interstitial complexes, is lower in samples without visible {311} de- 
fects. Interstitial clusters can become important sources for boron TED at low temperature 
at which the extended {311} defects are still stable against dissociation. Several defect 
states have been associated with interstitial clusters One of them is the P6 center 

commonly observed by electron paramagnetic resonance (EPR) measurements in low-energy 
ion-implanted, proton- or neutron-irradiated silicon [^Q. The P6 center has {100} symme- 
try, distinct from the typical {110} symmetry of vacancy- related defects. At low temperature 
(200 K), the symmetry of the P6 center is either C2 or Cih with a twofold axis parallel to 
the (100) direction. A thermally-activated symmetry transition occurs at room temperature 
(300 K) and the symmetry of the P6 center becomes 0- Motional averaging effects 
have been suggested to cause the transition 0,0. Complementary measurements of the 
^^Si hyperfine structure and the stress alignment indicate that the P6 center arises from 
di-interstitial defects. The stable di-intersitital is an important "precursor" of interstitial 
clusters and more extended defects such as {311} defects. Yet, the structure and dynamics 
of the di-interstitial have not been fully understood at an atomic scale and no first-principle 
calculations are available to our knowledge. 

In this letter, we propose a microscopic structure of a di-interstitial whose low- 
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temperature symmetry properties and electronic structure are consistent with the P6 center. 
Furthermore, the structure and transition paths between different defect orientations (Fig. 1) 
can account for the experimental symmetry transition to D2d at room temperature . The 
activation energy of the defect reorientation in our model is 0.5 eV in excellent agreement 
with 0.6 eV, the experimental activation energy extracted from the response of the P6 center 
to an external uniaxial stress [H. The experimental characterization of the donor level is 
also consistent with the calculated defect gap states. 

The atomic and electronic structure of di-interstitial defects are determined by ab ini- 
tio total energy calculations within the local density approximation (LDA) and generalized 
gradient approximation(GGA) The local minimum and metastable structures are 

obtained by structural minimizations using the conjugate gradient method. Starting geome- 
tries are stable configurations obtained by molecular dynamics simulations and structural 
relaxations using tight-binding Hamiltonians |[T0|-[T2|. In the ab initio calculations, the struc- 
tural minimizations employ two supercell sizes: (i) supercell A consisting of 120 atoms and 
(ii) supercell B consisting of 72 atoms. The total energy and structure are well converged 
with the plane- wave energy cutoff of 140 eV and 4 A;-points in supercell A [|T^. 

Figure 1(a) shows our di-interstitial model of Cih symmetry, with the twofold axis parallel 
to the z axis 0. The basic constituents of the model are a center atom Iq and dumbbell 
atoms /i-/2. Lee recently proposed a di-interstitial model composed of the same building 
blocks p|,p!^. Our model differs from Lee's model in (i) the orientation of the dumbbell 



atoms and (ii) the location of the center atom. The dumbbell structure, aligned parallel 
to the [110] direction, resembles the (110) interstitialcy, the most stable point defect [jl5|. 



Our suggestion that the di-interstitial can be formed when an interstitial is captured by the 
(110) interstitialcy is supported by the positive binding energy of the Cih di-interstitial with 
respect to isolated interstitials (Table I). 

Four equivalent Cih di-interstitial configurations, distinguished by the location of the 
center atom, can be constructed with three interstitials sharing one regular lattice site de- 
noted as in Fig. 2(a). The low temperature symmetry containing a twofold symmetry axis 



parallel to (100) is consistent with the Cih symmetry of our model @J^. Figure 2(b) illus- 
trates the four equivalent sites under symmetry which the center atom can occupy [|^. 
Thermal averaging motions among four Cih configurations with a small activation energy 
barrier can result in the room temperature D2d symmetry. 

Indeed, we can identify transition paths that lead to the thermal averaging motion be- 
tween four local minima. Three transitions - denoted as Tj, i = x,y, z, according to the 
twofold symmetry axis of its saddle point - have the same energy barrier of 0.5 eV. 
Tz transition : The transition from the center atomic site Jq to its mirror image 0, of, 
Fig. 1(a), with respect to the {110} plane can occur by a displacement of the atom Jq 
along the [110] direction. The saddle point of the transition has symmetry with a 
symmetry axis along z. The orientation of the dumbbell atoms remains the same after the 
Tz transition. 

Tx transition (Fig. 1(b)): A displacement of /i along the [101] direction results in a lo-h 
dumbbell pair and the atom I2 becomes the center atom denoted as Jq. The symmetry axis 
of the saddle point {C2V in Fig. 1(b)) is parallel to the [100] direction. The Ty transition 
is similar to the T^, transition, I2 moving along the [Oil] direction. As the result of the 
and Ty transitions, the orientation of the dumbbell atoms changes from the [110] direction 
to the [110] direction. 

The room temperature D2d symmetry can be explained by the thermal averaging mo- 
tion of three interstitials /o,/i and I2 alternatively occupying the four center atomic sites 
(Fig. 2(b)). Note that the Tj transitions leading to the thermal motional averaging of the 
D2d symmetry involve a displacement of only one of three interstitials along a C2V saddle 
point at each transition. The schematics of the potential surface along the high symmetry 
path is presented in Fig. 3. No other local minima are found along the path Cih ~ C2V — Cih- 
The experimental D2d symmetry is very plausible within our model, since the low transition 
energy barrier of 0.5 eV permits frequent T^, Ty and Tz transitions to motionally average the 
four Cih configurations accessible to three interstitials. 
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Structure and energetics of di-interstitial. Table I shows the formation energies Ef and 
binding energies Ei, of the Cih di-interstitial, the metastable di-interstitial and the (110) 
interstitialcy. Our LDA formation energy of the (110) interstitialcy agrees with that obtained 
by a previous LDA calculation ||T5[. In general, the bond lengths involved with the defect 
core atoms /o,/i and I2 are slightly longer than the bulk bond length of 2.35A. Charge 
density analyses confirm that the defect bonds are indeed weaker than the ideal bulk bonds. 
The strongest bond for both Cih and C2V di-interstitials is formed between the dumbbell 
atoms Ji and l2- 

Supercell B composed of 72 atoms is not large enough to quantitatively describe the 
energetics and structural properties of the Cih and C2V di-interstitials. We find a reduction 
in the transition energy barrier A from 0.6 (supercell B) to 0.5 eV (supercell A) due to the 
relaxation of the third- and fourth-neighbor-shell atoms from the defect core. The bond 
lengths associated with the defect core can differ by as much as 0.3 A between supercells A 
and B. The structural properties of supercell B are qualitatively different fTom those obtained 
by supercell A |T^. However, supercell A is sufficiently large for the converged energetics 
and structural properties |]TB|,|T^. 

We find large discrepancies between LDA and GGA in the formation Ef and binding 
energies Eb, obtained by supercell calculations with different number of atoms. This is at- 
tributed to the larger bulk cohesive energy of GGA compared to LDA. On the other hand, 
the relative stability of the Cih and C2V di-interstitials indicated by A and structural prop- 
erties are insensitive to the choice of the exchange-correlation potential V^c- The formation 
energy differences between supercells A and B are of the order of 0.1 eV within LDA. Simi- 
lar differences are obtained with GGA. The dynamical process for the symmetry transition 
schematically shown in Fig. 3 is not influenced by V^c- Typical differences between LDA 
and GGA results are less than 0.005 A in bond lengths and less than 0.5° in bond angles for 
the same supercell size. 

Electronic structure of di-interstitial. The experimental analysis of the EPR signal of the 
positively charged P6 center indicates that the donor level is strongly localized on the center 
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atom 0. Figure 4 shows that the defect gap states of the Cih di-interstitial are strongly 
locahzed. The donor level has strong jo-character mostly localized on the atom Jq, while the 
acceptor level is mainly localized on the dumbbell atoms /1-/2. The donor level of the Cih 
di-interstitial is located at Ey + 0.1 eV and the acceptor level at Ec — 0.2 eV. 

For the metastable C2V di-interstitial, the donor and acceptor levels become almost de- 
generate and form a deep level at + 0.4 eV. The metastable C2V di-interstitial becomes a 
stable defect when it is positively charged. The dashed line in Fig. 3 presents a schematic 
potential surface along the high symmetry path. Self-consistent calculations give a smaller 
energy difference of 0.14 eV between C^f^ and configurations. Our calculations on 
neutral and doubly-charged di-interstitials further support the effective symmetry of 
the P6 center even at room temperature. 

Link to extended defects. The stable di-interstitial is the link between point defects and 
extended defects. The dumbbell atoms of the di-interstitial are reminiscent of the (110) 
interstitialcy, suggesting the di-interstitial is formed by an existing (110) interstitial cy cap- 
turing an interstitial. Previously, we found a stability hierarchy of intestitial defects from 
molecular dynamics simulations with a tight-binding Hamiltonian |TD[: the formation en- 



ergy decreases in the order of interstitial clusters — >■ interstitial chains {311} defects. A 
strong tendency of clustering of interstitials has been also predicted from classical molecular 



dynamics simulations [0. When the interstitials are saturated such as in ion- implanted 
samples, the di-interstitials can play the role of nuclei for the extended interstitial defects, 
initially by providing interstitial sinks to form elongated interstitial clusters and eventually 
interstitial chains. We note that di-interstitials, larger interstitial clusters, and interstitial 
chains are composed of a common building block, the (110) interstitialcy |]10|. The aggre- 
gation of interstitial chains, in turn, can lead to the formation of extended {311} defects. 

In conclusion, we present first-principle calculations for structural and electronic prop- 
erties of di-interstitial defects and relate our stable di-interstitial structure to the P6 cen- 
ter. The Cih di-interstitial model and the transition paths involving the metastable C2V 
di-interstitial can account for different symmetries of the P6 center both at low and room 



temperature. The activation energy of the defect reorientation is 0.5 eV in excellent agree- 
ment with experiments. The localization of the donor level also agrees with the experimental 
characterization. 
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TABLES 

TABLE I. Di-interstitial {DI ) formation energies Ef and binding energies Ei,. The binding 
energy is defined as Ei, = —{Ef — 2Ef{I)), where Ef(I) is the formation energy of an isolated (110) 
/nterstitialcy. Supercell A consists of 120 bulk atoms, while supercell B consists of 72 bulk atoms. 
The orientation of the super cells is [110] x [110] x [001]. The energy barrier A is the total energy 
difference between Cih and C2v di-interstitials for the same supercell. 
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FIGURES 

FIG. 1. (a) Atomic structure of a Cih di-interstitial projected on the {110} and {110} planes 
and (b) atomic structures of the di-intcrstitial core during a transition between two Cih configu- 
rations. A regular lattice site located at the origin of the axes is shared by three atoms - (i) the 
center atom Iq and (ii) the (110) dumbbell atoms /1-/2. The site is related to the site Iq via a 
mirror reflection across the {110} plane. Three transition paths T^, Ty and Tz with the same energy 
barrier of 0.5 eV are identified. The C2v symmetry of the saddle point has a twofold axis parallel 
to the subscript of the transition Ti^i = x,y,z. For example, the transition in (b) displaces the 
center atom Iq to Iq site and the twofold axis of the saddle point in the middle is parallel to the 
[100] direction. 

FIG. 2. (a) Atomic structure of a Cih interstitial and (b) schematics representing the four 
equivalent sites of the center atom Iq under D2d symmetry operations. The regular site indicated 
by o is shared by three intcrstitials Iq, Ii and l2- The atomic indices are assigned based on the 
distance from the origin and the symmetry to aid visualization. The experimental symmetry 
of the di-interstitial at room temperature can be explained by the thermal averaging motion of 
Iq,Ii and I2 alternatively occupying four equivalent sites. As the result of a Tj transition, the 
center atom Iq "effectively" moves to the corresponding site as indicated by the arrows. Note that 
black atoms lie below the grey atoms along z axis. 

FIG. 3. Schematic potential surfaces along the highly symmetric path of Cih — C2V — Cih 
transitions for the neutral state (solid line) and for the positively charged state (dashed lines) at a 
given chemical potential /x. The neutral C2V di-interstitial is metastable and the energy barrier A 
is 0.5 eV. In contrast, the positively charged C2V di-interstitial is stable and the energy difference 
5 between Cih and C2v di-interstitials is only 0.14 eV. 
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FIG. 4. (a) Isosurfaces of gap state amplitudes of the Cih di-interstitial, (b) local density of 
states (LDOS) on the center atom Iq (solid line) and the dumbbell atoms I1-I2 (dashed line) of 
the Cih di-interstitial and (c) those of the C2V di-interstitial. The reference energy of (b) and 
(c) is Ey, the valence band maximum of bulk silicon in the 120-atom supercell A. The calculated 
bulk conduction band minimum, Ec, is indicated by an arrow and is about 0.8 eV from E^. The 
energy gap of the Cih di-interstitial is about 0.5 eV, while the donor and acceptor levels of the C2V 
di-interstitial are almost degenerate. LDOS of the gap states is highlighted by the black (gray) 
areas for the donor (acceptor) level. 



12 




Kim et al.^ Fig. 1 




Kim et a/., Fig. 2 




Kim et al.^ Fig. 3 




Kim et al.^ Fig. 4 



